Fixing Two-Nucleon Weak-Axial Coupling L^a From /i d Capture 
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We calculate the muon capture rate on the deuteron to next-to-next-to-leading order in the 
pionless effective field theory. The result can be used to constrain the two-nucleon isovector axial 
coupling Li,a to ±2 fm 3 if the muon capture rate is measured to 2% level. From this, one can 
determine the neutrino-deuteron break up reactions and the pp-fusion cross section in the sun to a 
same level of accuracy. 
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The strong evidence of neutrino oscillations observed 
at the Sudbury Neutrino Observatory (SNO) l] is based 
on detecting the 8 B solar neutrino flux through the fol- 
lowing three reactions: 



v e + d^p + p + e~ (CC), 
v x + d^p + n + v x (NC), 
v x + e~~ — > v x + e~ (ES). 



(1) 



The charged current (CC) reaction involves only the elec- 
tron neutrinos, while the neutral current (NC) reaction 
and elastic scattering (ES) involve all the active neutri- 
nos (x — e,n,r). The v e and v x fluxes are found to be 
significantly different |jj . Further detailed measurements 
of the fluxes could sharpen the constraints to neutrino 
oscillation parameters and provide precision tests to the 
standard solar model 0. However, while the ES cross 
section is known to high accuracy, the CC and NC cross 
sections have hadronic uncertainties. As shown by But- 
ler, Chen and Kong 3], the dominant uncertainties in 
low energy CC and NC cross sections comes from the cou- 
pling of a two-body isovector axial current, L±,a, in pion- 
less effective field theory (EFT(yf)). The potential model 
results of Refs. and can be reproduced by different 
choices of L\,Ai indicating that the ~ 5% difference be- 
tween the models comes from the different assumptions 
about the short-distance nuclear physics. There are other 
interesting weak reactions involving the same two-body 
current, for example, the pp and pep fusion processes 
(pp — > de + v e , ppe~ — ► dv e ) which power the sun 
It is one of the great current interests to measure these 
neutrino fluxes to further test the standard solar model. 

Recently much effort has been going into determin- 
ing the effective two-body axial current interaction 
Butler, Chen, and Vogel attempted to fix L\ : a from re- 
actor antineutrino-deuteron breakup reactions, and they 
found L 1 a — 3.6 ± 5.5 fm 3 ||. Chen, Heeger, and 
Robertson obtained L x .a = 4.0 ±6.3 fm 3 by using SNO's 
CC and NC data, calibrated by the ES events of SNO 
and Super-Kamiokande(SK) [9|. Schiavilla et al.'s idea 
[loj ] of using the tritium (3 decay rate to control the 
strength of the two-body current was adopted by Park 
et al. in their hybrid EFT calculation [llj, and the pp 



fusion rate was predicted with a small error. When com- 
pared with the KFT(f) calculation 0, their result yields 
L x ,a = 4.2 ±2.5 fm 3 . 

In this paper we aim to make a high-precision deter- 
mination of Li^a from the \x~ d capture process 



(2) 



by calculating the rate to next-to-next-to leading order 
(N 2 LO) in EFT(yf). The fi~d capture rate has been mea- 
sured previously by different groups with rather different 
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results r cx P = 470 ± 29 s" 1 12] and T cx p = 409 ± 40 s" 1 
[l3j . A measurement of this rate with 1% precision 
is under investigation at PSI 0]. An earlier poten- 
tial model calculation [l5| gave T = 397 ~ 400 s _1 . 
More recently, the h ybr id approach mentioned above 
gave r = 386 ±5 s" 1 [l|. 

A concern in applying EFT(/r) to the d capture is 
that the energy transfer into the hadronic system might 
be too large to apply EFT(/r). However, as shown in 
Ref.0 and also in this calculation, the contribution to 
the total rate from high-energy neutrons is small, and 
it is possible to impose a neutron energy cut to isolate 
the low-energy (< 10 — 20 MeV) neutron events without 
significantly increasing the statistical errors 0, ^| . 

Effective field theory is useful when low and high en- 
ergy scales in the problem are widely separated. For low- 
energy processes, short-distance physics can be taken into 
account by local operators in an effective lagrangian in- 
volving only low-energy degrees of freedom. For v(v) — d 
scattering with neutrino energy below 20 MeV and /i~d 
capture with small final-state neutron energy, the pion 
and other meson exchanges are not dynamical degrees of 
freedom, and their physics can be captured by contact in- 
teractions involving nucleons and the external currents. 
To make predictions with controlled precision, calcula- 
tions are done with the perturbative expansion param- 
( 1 s ) 

eter Q = (l/<4«°',7 , p) / A , which is the ratio of light 
to heavy scales. The light scales include the inverse S- 
wave neutron-neutron scattering length l/a nn ° = —10.6 
MeV in the channel, the deuteron binding momen- 
tum 7 = 45.7 MeV in the Z S\ channel, and typical nu- 
cleon momentum p in the system. The heavy scale A is 
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set by the pion mass m^. This EFT (ft) (see e.g. \L7\) 
and its dibaryon version 0, 0, [2(| have been applied 
successfully to many processes involving the deuteron, 
including electro-ma gnet ic p rocesses such as Compton 
scattering 7c? — > 7c? |2lL |22|. n» — * c^y relevant to the 
big-bang nucleosynthesis [23|, 123. weak processes such 
as vd reactions for SNO physics Q, the solar yy fusion 
process 0, 0j , and parity violating observables . For 
reviews on three-body systems, see |2fij |. 

The effective Lagrangian for the CC weak interaction 
is given by C cc = —^uJiGpVudl+J^/VZ + h.c, where 
Gf = 1-166 x 10 _5 GeV -2 is the Fermi coupling constant 
and u>i = 1.024 takes into account the inner electroweak 
radiative correction j^. = F^7^(l — 75)^ is the lep- 
tonic current. The quark current J~ = V7 — A~ — 
(V} — A 1 ) — i(V? — A 2 ^ ) contains both vector and axial- 
vector interactions, where the superscripts 1 and 2 are the 
isospin indices. At the scale relevant to nuclear physics, 
the quark current need be matched to a hadronic current 
which in general contains one-nucleon, two-nucleon, etc., 
operators. 

Up to the order of our interest, the one-nucleon isovec- 
tor vector and axial vector currents are 
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where V = V + V , d = d a + d , and V = V - V . The 
superscript a is the isospin index, fi^ = (fj, p — /i„)/2 = 
2.353, is the isovector magnetic moment. Isovector Dirac 



charge radius (r 



2\i = l 



,2\P 



2\I=1 
A 



0.873 fm 2 , and 
~ 0.45 fm 2 . We 



the isovector axial-charge radius (r 

have neglected terms of order p 2 /Mjy or even n^p 2 /M^. 
The pseudoscalar form factor is, to a good approxima- 
tion, dominated by the pion-pole G p (q 2 ) — ^pz^p whose 

q 2 dependence will not be expanded because the momen- 
tum transfer |q| is of order muon mass M M with low en- 
ergy final state neutrons. The G p contribution to the 
axial current is counted of order Q. 

The lowest dimensional two-nucleon isovector currents, 
in the dibaryon version of EFT (ft), relevant to the fi~ d 



capture process are 



V (2) 



,11, 



ekijt\VjS a + h.c. , 
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N 
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db 



(t{s a (4) 

where ti and s a are dibaryon fields for the two-nucleon 
3 Si and 1 So states, respectively. The second term in 
A 



^2" is induced by the G p term in the one-nucleon cur- 



( 1 s ) 

rent. r { ° lj = 1.764 fm and r„„ = 2.8 fm are the 
effective ranges in triplet and two-neutron singlet chan- 
nels, respectively. Lf b and Lf b A are coupling constants 
in dibaryon formalism. The vector current is N 2 LO and 
its coupling Lf b = —4.08 fm has been determined by the 
rate of n + p — > d + 7 near threshold. The axial cur- 
rent is NLO, and its coupling Lf A is proportional to the 
renormalization-scale-/i-independent Li t A in Ref. [jj as 

L \,A 



^-Li_a, through which Lf b A is related to the /it- 
dependent Li^A(pf) in Ref- lil- The numerical relation be- 
tween Lf A and Li tA (n) is Lf A = -13.8 + 0.28Li,^(M 7r ), 
where L± t A(M w ) is in units of fm 3 and has a natural size 
- 6 fm 3 . ' 

The [i~d atom has a ground state with a hyperfinc 
structure, corresponding to the total angular momentum 
F = 3/2 and F = 1/2. The /j,~d capture process is 
known to take place almost uniquely from the doublet 
F = 1/2 state. The differential capture rate for muon 
and deuteron in their specific polarization states can be 
written in terms of leptonic tensor and hadronic ten- 
sor W av as 



dT(S,j) = u t G 2 F \V ud \ 2 E\iP(0)\ 2 
dEdn 32^ 2 M M (1 + f^) 



l^(S)W^), (5) 



where the 
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em Mu+Md 



is the lS-state 



wave-function-at-origin-squared, E is the outgoing neu- 
trino energy, and Md is the mass of the deuteron. 

The capture rate depends on the polarization vector of 
the muon S M and the deuteron polarization vector £. The 
leptonic tensor is given by P"(S) = A(k^k' v + k'^k v - k ■ 
tfgnu + ieWkpk'v) - 4M AI (S^/c" / + k'»S v - k' ■ Sg^ v + 
ie^P a S p k' a ), where k = (M M ,5) and k' = (E,k'), with 
E = \k'\, are the four- momenta of initial muon and final 



respectively. The hadronic tensor is 



W^ v (i) = -lm 

TT 



d 4 xe lqx T(d\J+(x)J-(0)\d) 



(6) 



where = k^ — k', and \d) = \d(P, £)) is the deuteron 
state with momentum P = (Md,0) and polarization £. 
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The diagrams contributing to W^^) up to N 2 LO are 
shown in Fig. 1. A straightforward calculation finds 



dT_ = G 2 F \V ud \ 2 E 2 \^(0)\ 2 
dE 
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(7) 



-2 X /„,2\J=X„2 
6 ' 

-2 1 /^2\/=l„2\ 



1 + e z ^(r' z ) L ~ Y q 2 is the Dirac form factor, 



where Gy 

and = <7a(1 + £ 2 g( 7 " 2 )^ =1 <? 2 ) is the axial form factor. 
The e is formally introduced to keep track of the Q ex- 
pansion. After expanded in e and truncated at e 2 , the 
N 2 LO result is obtained by setting e = 1. The functions 
F la ,x6, F 2 , and F 3Qi3bi 3 C are from diagrams (a), (6), and 
(c), respectively, in Fig. 1 
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with the re-scattering amplitude in the singlet channel 
as A( ls °Hp) = -^(-^T " 4^ S< V + ip)~\ and 
functions G12 are defined as Gi = tan -1 ( 2 / T ^ ip \ ) + 



-Lf b A and G 2 



tan x (t 



The 



energy injection into the two-neutron system is = 
Mfj, — E — (M n — M p ). The relative momentum be- 
tween the two final-state neutrons is 2p, with p = 
^/MnAe-J 2 - E 2 /A. 

Power counting for the present calculation is rather 
tricky, and it would be misleading, for example, to just 
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FIG. 1: Diagrams contributing to U P to N 2 LO. The 

lines inside the loop are the nucleon propagators, and the 
thin and thick lines outside represent the triplet and singlet 
dibaryon fields, respectively. The small open circle represents 
the insertion of hadronic one-body current. In diagram (c), 
the large gray circle indicates the two possible hadronic cur- 
rent insertions shown in the second row where the dark square 
denotes the two-body current associated with Li,a and L\. 




E nn {MeV) 

FIG. 2: (Color online) The differential capture rate dT/dE„ n 
calculated using = 6 fm 3 . The doted line is the LO re- 
sult. The dashed and solid lines, which sit on top of each other 
in this scale, are the NLO and N 2 LO results, respectively. 



use the p/m^ to estimate the accuracy of the expansion. 
For example, it is well-know that in the nucleon-nucleon 
scattering, the effective theory without pion works rather 
well at the nucleon momentum on the order of 100 MeV, 
a value close to the pion mass. This can also be seen in 
the present calculation because the NLO result at large 
neutron momentum does not completely modify the lead- 
ing order result, whereas the naive power counting would 
indicate otherwise. Fortunately, for muon capture, the 
most of the events occur when the neutron momentum is 
about half of the muon mass, that is, about 50 MeV. 

In Fig. 2 we show the differential rate dT/dE nn 
in terms of the relative motion energy E nn = 
2{^/M N + p 2 — Mn) of two final-state neutrons in the 
region where the EFT(/f) calculation is most reliable. It 
is clear from the figure that the differential rate in the 
energy region E nn > 10 MeV is very small, and is neg- 
ligible for E nn > 15 MeV. By comparing the results of 
LO, NLO, and N 2 LO, we find good convergence of the 
expansion. 

In the case that a neutron energy cut can be imposed 
on experimental data Q|, it is possible to define and 
measure the integrated capture rate up to a threshold 
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FIG. 3: (Color online) The integrated capture rate V(E, 
calculated using L\,a = 6 fm 3 for the relative energy E^ n of 
the two neutrons up to 30 MeV. The description of the lines 
is the same as that in the caption of Fig. 2. 
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(9) 



The result up to 30 MeV is shown in Fig. 3. In the 
whole energy region, the NLO contribution is less than 
10% of the LO contribution, while the N 2 LO contribution 
is less than 1%. This small size of N 2 LO contribution 
is accidental and does not happen for unpolarized and 
F = 3/2 rates. For example, the unpolarized rate has 
the expansion (for Lx,a = 6 fm 3 ) 



range of £™ , the size of b(E™) is about 1.3 - 1.5% of 
the size of a(E r ^ T \). This shows how an error in capture 
rate is translated into an uncertainty of £i,A- 

Table 1: 

Coefficients functions a(£^) and b(E™ l ) for specific val- 
ues of two- neutron relative energy E^ n from the EFT(7f) 
calculation. 



£&(MeV) 


5.0 


10.0 


15.0 


20.0 


«(K h n)(s- 1 ) 


239.2 


308.0 


332.0 


342.3 


fe^Hs-ifm- 3 ) 


3.3 


4.2 


4.7 


4.9 



In summary, we calculated the fi d capture rate us- 
ing EFT(tt) to N 2 LO. The major goal is to fix the two- 
nucleon isovector axial coupling constant Li a from fu- 
ture precision experimental data. An experimental result 
on the integrated rate up to some neutron energy E^ 
with a 2% error should be able to, through comparison 
with our calculation with theoretical error 2-3%, fix the 
L\,a with error ±2.0 fm 3 . This in turn allows us to de- 
termine the neutrino deuteron breakup cross section and 
the pp fusion rate in the sun to 2-3%. 

This work was supported by the U. S. Department of 
Energy via grant DE-FG02-93ER-40762 and the NSC of 
Taiwan. We thank S. Ando, D.W. Hertzog, P. Kammel, 
K. Kubodera, and T.-S. Park for helpful discussions. 



- unpol. 



1^.(1 + 5.3% + 4.9%), 



(10) 



where the NLO correction is abnormally small, and 
NNLO is of the normal size. A similar expansion is ob- 
tained for Li.a = —6 fm 3 : 



unpol. 



v LO 
unpol. \ 



10.9% + 5.2%) 



(11) 



which shows a nice convergence pattern. Based on this 
trend, we assign a 2-3% correction at NNNLO, corre- 
sponding to an error in L\ a ~ 2 fm 3 . This is consistent 
with the naive estimation of 3% if the small expansion 
parameter is 1/3. A calculation shows the N 3 LO final 
state P-wave re-scattering contributes only ~ 1%. Fur- 
thermore, the result is insensitive to the uncertainty in 

f ^ S ) 3 

a-nn ■ Choosing Li,a = 5.6 fm , the energy dependence 
of our result matches the previous hybrid calculation very 
well 0. 

To extract Li,a from experimental data, it is useful to 
provide the dependence of the rate on L\ t A 
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HEn h n)Ll,A, 
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where Li,a is in unit of fm 

nth \ _„ J U T^th 



The energy-cut dependent 

nth , 



functions a(E™ n ) and b{E™ n ) for a set of E™s are listed 
in the Table 1, from which we observe that, for the whole 
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